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ABSTRACT 

In this paper, we present a redshift space reconstruction scheme which 
is analogous to and extends the Perturbative Least Action (PLA) method 
described by Goldberg & Spergel (2000). We first show that this scheme is 
effective in reconstructing even nonlinear observations. We then suggest that by 
varying the cosmology to minimize the quadrupole moment of a reconstructed 
density field, it may be possible to lower the errorbars on the redshift distortion 
parameter, f3 as well as to break the degeneracy between the linear bias 
parameter, b, and Qm- Finally, we discuss how PLA might be applied to realistic 
redshift surveys. 

1. Introduction 

The current generation of galaxy redshift surveys is producing an avalanche of data 
about the structure of the universe. The Sloan Digital Survey (SDSS; York et al. 2000) has 
already measured the redshifts of ~ 30, 000 galaxies (X. Fan, private communication) of the 
~ 10 6 galaxy redshifts over 10,000 deg 2 which the survey will ultimately cover. The Two 
Degree Field redshift survey (2dF; Colless 1999) has measured redshifts for ~ 10 5 galaxies, 
and will eventually measure a quarter of a million galaxies over 2000 deg 2 in the southern 
hemisphere. Meanwhile, the IRAS 0.6 Jy Point Source Catalog redshift Survey (PSCz; 
Saunders et al. 2000), with about 15,000 galaxies over virtually all of the sky, provides a 
fertile testbed for cosmological models and methods. 

As impressive as these surveys are, they are limited to providing a somewhat distorted 
snapshot of the universe. For example, it is both suggested observationally (Hubble 1936; 
Oemler 1974; Davis & Geller 1976; Kaiser 1984; Santiago & Strauss 1992; Blanton 2000 
and references therein) and predicted (Davis et al. 1985; Bardeen et al. 1986; Blanton 1999 
and references therein; Dekel & Lahav 1999) that the luminous structure, that which the 
surveys record, is biased with respect to the underlying matter density of the universe. 
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Additionally, though the Hubble relation can be used to give an approximate 3-dimensional 
picture of structure, peculiar velocities (e.g. Strauss & Willick 1995 and references therein 
for a review) cause a distortion of the structure along the line of sight. Since the peculiar 
velocity field is a function of the underlying mass field and cosmology, velocity distortions 
and bias are intimately related. 

Perturbative Least Action (PLA; Goldberg & Spergel 2000, hereafter GS) was shown 
to be an excellent technique for the reconstruction of nonlinear structure in real space. In 
this paper, we extend PLA into redshift space, and show that it is an excellent tool for 
extracting information from redshift surveys even into the nonlinear regime. 

However, before going too far afield, it will be useful to review some of the basic 
issues involved in redshift space distortions of the density field, and define some of the 
symbols which will be used throughout this paper. This discussion is not meant to be 
comprehensive, however, and the interested reader will certainly benefit from some of the 
excellent reviews on the subject (Hamilton 1998; Hatton & Cole, 1998; Zaroubi & Hoffman 
1996; Strauss & Willick 1995; Sahni & Coles 1995; Kaiser 1987). 



1.1. Definitions and Conventions 

Let us consider an observer in an expanding universe. Hubble's law states that in a 
Friedman-Robertson- Walker universe, the distance, d to a test particle with redshift, z, will 
be, to first order in z: 

cz ~ H d FRW , (1) 

where H is the Hubble constant at the present epoch, d is the distance to the particle, 
and the approximation comes from the fact that relativistic effects become important at 
high redshifts. However, we will confine our discussion to the non-relativistic case and this 
approximation throughout this discussion. 

Hubble's Law assumes that a particle (galaxy) is at rest in comoving coordinates. A 
particle with a local, peculiar velocity, v, will have a redshift of: 

cz = H d + x • v , (2) 

where x points along the line of sight of the test particle. Since it is actually this redshift 
that we observe, and not the position of the particle, it is worthwhile to construct a 
comoving coordinate which reflects the observations of the observer at the origin. We define 
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the comoving redshift space coordinate: 

( ax ■ x\ 

where we have substituted ax for v and ax for dx in our implicit definition of the peculiar 
velocity. Comparison with equation (H) yields the relation, s(t) = cz/a(t). 

Equation (|3|) shows that the mapping of x to s is inherently non-invertible. Any 
trajectory of x will yield a single trajectory in redshift space, but the converse does not 
necessarily hold. 

Indeed, even in the Zel'dovich regime (defined below), the infall of matter from both 
sides (front and back) of a structure can give rise to a "triple-valued zone" (see Strauss 
& Willick 1995, §5 for a discussion), as matter at different physical distances from the 
observer appear to have identical redshifts due to conspiracy between the Hubble flow and 
the peculiar velocity. We will attempt to disentangle this degeneracy in § 0. 



1.2. The Linear Regime 

For now, let us consider the linear regime. By convention, particle i sits at position, qj, 
at a = 0, and that the ensemble of q« forms a uniform grid. If a field is linear, that is, if 

'<».«) w 

remains small at all times, and if the velocity field initially has no curl, then the Zel'dovich 
approximation (Zel'dovich 1970) can be used to give the trajectory of a particle as: 

Xi(t) = qi + D(t) Pi , (5) 

where D(t) is a cosmology dependent, monotonically increasing growth factor, normalized 
to unity at the present, and p» is the final displacement of particle, i, from its initial 
position. 

Fields for which equation (HD well approximates the trajectories of all particles at all 
times will be referred to as Zel'dovich fields. This is to be distinguished from fields in the 
"linear regime" for which: 

5(qi,t) = -D(t)V q -p(qi) (6) 

holds at all times. When perturbations are very small, both equalities hold. However, as 
perturbations become larger, equation (HD breaks down first, and densities evolve according 
to a more complex function of time. 
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In the Zel'dovich regime, the redshift space coordinate evolves as: 



s(t) {zel) = q + D(t)p + x 



(a(t)D(t) N 



where we define f(&M, &A,t) such that: 



q + D(t)p + x [D(t)f(Q M , fi A , t)x • p] (7) 



a(t)D(t) 
a(t)D(t) 



For = 1, this function is a constant in time. At t = to? a good analytic fit can be given 
by f(Q M ,Q^,t ) ~ f2^ 6 (Peebles 1980). This function is normally used in discussions of 
bias, and is generally combined with the linear bias parameter, b, to relate the divergence of 
the velocity field to the overdensity of the galaxy field at the present day, via the parameter, 
(3 = /(VIm, &a, to)/b, where 5 g = b5 for a linear biasing model. 

In this section, we will be treating only unbiased fields, and have introduced j3 as a 
means of comparing this discussion of redshift space distortions to the standard approach 
(e.g. Strauss & Willick 1995). In §3, we'll return to the degeneracy in (3 and show how PLA 
might be used to break it. 



1.3. The Distant Observer Approximation 

Up to this point, we've treated linear redshift space distortions with more or less full 
generality. However, since our ultimate goal is to apply these distortions in the context of 
PLA, we will want to make some simplifying assumptions. For example, PLA (GS) uses 
a Particle Mesh (PM) Poisson solver (Hockney & Eastwood 1981). This method takes 
advantage of Fast Fourier Transforms (FFTs), which assume Cartesian coordinates. 

The general form of the comoving redshift coordinate, s, above (equation |3[), is 
not separable in Cartesian coordinates. If this form were to be applied in general, one 
would wish to describe coordinates with spherical harmonics (e.g. Susperregi 2000). 
The primary purpose of the current discussion, however, is to examine the underlying 
dynamics, and while treating redshift space distortions of nearby systems is undoubtedly 
of cosmo logical interest, the matter at hand is greatly simplified by assuming the distant 
observer approximation (d.o.a.). 

In the d.o.a. we essentially assume that the system of interest is sufficiently far away 
that the x, are parallel for all particles. For convenience, we will label the third orthonormal 
coordinate (the z-axis), as the line of sight. Using this definition, we redefine the comoving 
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redshift coordinate: 

s a = x° + 5«^ , (9) 
a 

where a denotes the index of the direction vector, and 5^8 is the Kronecker-5 function. 
In the linear regime, this becomes 

s a (t) = q a + D(t) P a [i + 6^f(n M ,n A ,t)] . (io) 

Thus, if a particle is observed at sf and its initial position, q i; and the cosmology are 
known, this expression may be inverted and combined with equation (|5|) to give: 

D(t)(s aF - a a ) 
i 1 + <5a 3 /(i2 M ,S2 A ,to)J 



Equations (lOjIl) may be combined to extend linear theory into redshift space. One 
must keep in mind that we have assumed that we know both the initial and final position of 
a particle in the comoving redshift space coordinate. As pointed out above, without both 
constraints, inverting the redshift coordinate becomes ill-posed. 



2. Method: Least Action in Redshift Space 

Even if we have the idealized set of observations discussed above, and have a complete, 
unbiased mapping of the density field in redshift space, as perturbations become large, 
complications will arise both in mapping the redshift space field to a real space one, and 
in reconstructing an initial density field. A number of researchers have attempted worked 
on the dual problems of reconstructing an underlying real space CDM density field from a 
velocity field, and the calculation of a velocity field field from a redshift survey. 

For the former problem, the POTENT algorithm (Bertschinger & Dekel 1989; Dekel, 
Bertschinger & Faber 1990; Dekel et al. 1999) uses the Zel'dovich approximation to relate 
a redshift/distance survey, in which one component of the peculiar velocity can be directly 
computed, to an underlying mass density field. Nusser et al. (1991) uses a nonlinear 
generalization to extend this into the nonlinear regime. Others (Kudlicki et al. 1999; 
Chodorowski et al. 1998; Chodorowski & Lokas 1997; Bernardeau 1992) use higher order 
perturbation theory to compute the relationship between the velocity divergence and real 
space density field. In particular, Chodorowski & Lokas (1997) point out that application 
of these methods may be used to break the degeneracy between the linear bias constant, b, 
and cosmology. 
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A related problem concerns the calculation of the peculiar velocity field from a galaxy 
redshift survey. Nusser & Davis (1994) use a quasi-linear correction to the Zel'dovich 
approximation to relate the redshift space density field to the peculiar velocity field. Like 
the methods listed above which use perturbation theory to related the velocity and density 
field, Chodorowski (2000) uses 2nd- and 3rd- order perturbation theory to expand and 
compare the real space and redshift space density fields. 

Others have taken a slightly different approach, which attempts to essentially solve 
these two problems simultaneously. Giavalisco et al. (1993) suggested that the Least Action 
approach described by Peebles (1989) could be used in redshift space with only a canonical 
transform of the coordinates. Schmoldt and Saha (1998) consider the difficulties of running 
least action reconstruction in redshift space, and test this by reconstructing the velocity 
field of the Local Group. 

Susperregi & Binney (1994) modified this approach and described a technique whereby 
one could expand the density and velocity fields in Fourier space, and write down the Least 
Action equations in Eulerian form. Susperregi (2000) took this a step further, and applied 
a similar code (using Spherical Harmonic transforms) to the IRAS 1.2 Jy redshift survey 
(Fisher et al. 1995). Each of these techniques use a smoothing filter on the density field 
such that they are in the mildly nonlinear regime at the present. 

PL A has several distinct advantages over these approaches. First, since they are 
inherently Eulerian, as perturbations becomes large, they no longer fairly sample the matter 
field. PLA, on the other hand, is Lagrangian in the sense that it performs the time integral 
over the particle field, rather than the density field. Moreover, the Eulerian PLA approaches 
assume a locally curl-free velocity field at all times, by the velocity-density relationship. 
While PLA generally assumes curl-free initial conditions, vorticity is permitted to develop. 

Ultimately, we want to reconstruct the underlying real space CDM density field and 
evolution from a set of observed galaxy redshifts. The approach taken in this paper differs 
from those discussed in GS in that we now deal with quasi-linear structure, rather than the 
highly nonlinear constraints. While previously, we were content to get a "realistic" set of 
initial conditions, here, our aim is reconstruct the details of the observations exactly. By 
doing this, we hope to disentangle degeneracies in bias, get a handle on the true cosmological 
power spectrum of perturbations, and examine the growth of large scale structure. 

In this section, we describe applying the PLA approach to redshift space constraints. 
In particular, we will deal with two main issues: 

I. Given some observed redshift space density field, S F (s), find final particle constraints, 
{sf }, which satisfy the density field, which maps to a corresponding initial uniform 
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field, {q^}, such that the constraints can be most easily satisfied. 

2. Given an initial particle constraints, {qi}, and final particle constraints in redshift 
space, {sf }, find the trajectories of particles which self-consistently satisfy the 
boundary conditions. 

2.1. Computing Particle Constraints 

We must first determine the boundary conditions for the particles' positions at a = 
and a = 1. Let's consider a set of idealized observations, in which a smoothed, complete, 
and unbiased density field in redshift space, S F (s), has been observed. 

In order to determine the final boundary condition, we generate a set of particle redshift 
positions which yield this observed density field. We begin by assuming by distributing 
particle positions uniformly distributed on a grid, Sj = q^, where the tilde over the redshift 
coordinate will be explained shortly. 

From here, we iterate in the following way. In each iteration, we take the density 
field of the current value of the particle positions yielding S(s). We then transform this 
density field into the target field by using the Jacobian determinant to create a laminar 
flow. That is, one can adjust the particle positions in the former grid by using a coordinate 
transformation: 

s^ = s, + V^(s). (12) 

In that case, the density field as measured in the primed frame compared to the unprimed 
frame will be: 

[1 + 5(s')] dV = [l + 5(s)]rf 3 s . (13) 

From the form of the transform, the determinant of the Jacobian is easily computed. For 
small perturbations: 

[l + 5(s')] ~ [l + 5( s )](l+V 2 ^) . (14) 

Assigning S(s') = S F (s), and using standard Fourier techniques, one can compute the 
scalar field, i/)(s), and taking the gradient, one can compute the coordinate transform, and 
therefore, can iteratively produce a particle map which satisfies the target density field. 

Since we assume throughout that structure evolves out of an initially uniform density 
field, the initial constraints on these particles must be the uniform grid, q^. This is 
exactly analogous to mapping from an unperturbed simulation since a uniform particle 
field run through an N-body will remain uniform, and hence serve as a perfectly legitimate 
unperturbed simulation. 



- 8- 



But which final positions correspond to the values of q^? In other words, given any set 
of particle redshift positions, sf which satisfy the density field, S F (s), what is the "best" 
permutation matrix, My, such that 



S; 



E M ii*I ? ( 15 ) 



We define Mjj as an N p x N p matrix which has exactly one "1" in each row and column, 
and zeros elsewhere. 

We need to define what we mean by "best". In general, we mean that we wish 
to compute boundary constraints which most naturally provide us with physically well- 
motivated orbits. Since the laminar flow method of generating final constraints necessarily 
assumes no shell crossings, for small perturbations, the permutation matrix will simply be 
the identity matrix. In this ultra-linear case, in which perturbations are assumed to be 
so small that there are no orbit crossings even in redshift space, we estimate the physical 
particle displacement from the uniform field as: 

v a = (16) 

The final physical position of particle, i, is thus related by equation (g). Taking the 
smoothed density field of Xj gives an estimate of the density field in real space. 

However, in many cases of interest, even where there are no shell crossings in real 
space, there are orbit crossings in redshift space. These are the famous "triple value zones," 
(Strauss & Willick 1995 §5.9) so named because a particle at a given redshift is a triply 
degenerate function of distance. Even if a system is dynamically in the linear regime and 
the real space density field is known with great accuracy, this degeneracy can arise. 

In order to illustrate this, we have run a simulation in which a plane wave density 
field is laid down along the line of sight. Its amplitude is such that there are triple- valued 
zones at a = 1. The solid lines in Figure [3] show the initial and final density field of this 
distribution in both real and redshift space. The solid line in Figure |2] demonstrates the 
existence of triple value zones. There, we plot the relationship between real and redshift 
particle coordinates for our simulated Zel'dovich pancake. 

Given the "observed" redshift space density field, 8 F (s) in panel d) of Fig. [I], we 
have already described how to generate a set of redshift space coordinates sf . Since the 
simulation is approximately in the Zel'dovich regime, it is assumed that if we can estimate 
the real space coordinates, we can use that information to estimate the position of a particle 
on the uniform grid,qj. However, given the degeneracy of the mapping from redshift to real 
space, this is no simple task. 
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P(xf|sf)cx[l + ^(x l )]^ 



X « X ( S i i 



Willick et al. (1997; also Sigad et al. 1998) suggest a likelihood approach to breaking 
this degeneracy, called VELMOD. Part of VELMOD relates the observed redshift space 
density field to a test value of the real space density field. We use a similar approach 
here. We first assume that the real space particle field is the one generated from the 
assumption of no orbit crossings. That is, equations fll6| ) is used to approximate the 
positions of the particles. One may then compute a real space density field from the particle 
field approximations, and from that the potential field, 0(x) can then be computed. The 
assumption of linearity gives the following relation: 

pW = -rifhn • < 17 > 

{a D + a D ) 

This makes s F a straightforward function of x F via equation (^j). For triply valued sf , we 
can then determine a posterior probability that the particle is at real space position, x F : 

OSi 

where 5 OLD (x) is the previous iteration of the estimated real space density field. This is 
simply the discrete form of the continuous distribution function used in the VELMOD 
approach. Using this distribution function, a real space coordinate is randomly assigned to 
each triple-valued particle. The real space density field and potential are then recomputed 
and the process is repeated until satisfactory convergence is reached. 

We apply this VELMOD-like approach to our Zel'dovich pancake. The dotted and 
short dashed lines in the lower-left panel of Figure [I] show the initial guess and final 
estimate of the real space density field. While the fit between the true and estimated real 
space density fields are good, they are not perfect. One way of thinking about this is that 
equation ( |17|) assumes that the final velocity of a particle is linearly proportional to the 
force on that particle. In the Zel'dovich approximation, this is a good assumption. However, 
as the limits of that approximation are approached, the relation between acceleration and 
velocity may become more complex. On scales on which shell crossings occur, the two may 
even be of opposite signs. One approach that people have historically used is to simply 
smooth density fields until all structures are linear. It should be noted this density field is 
not our final estimate of the real-space field. Rather, we are using it as a working model to 
set up redshift space constraints for PLA, which makes no assumptions about the linearity 
of orbits. 

This artificial steepening is also apparent in the estimated relationship between real 
and redshift space coordinates, as is illustrated by the open square points in Figure |[ 
Notice, however, that we have qualitatively reproduced the features of the tripled value 
zones. 
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At the final iteration of VELMOD, not only are assumed final positions of the particles 
computed, but so are their velocity/displacement vectors, p^. If linear theory approximately 
holds, then each of those particles ought to have originated at: 

q? = sf a - [i + 5*f(n M , n A , to)] Pt ■ (19) 

We may now return to the problem posed at the beginning of this section: How do we 
compute the "best" permutation matrix, My? We find the matrix which minimizes: 

X 2 \Mijqj, Qi] = E - MijQj] 2 • (20) 

i 

In order to actually perform this minimization, we apply a simulated annealing method 
(Press et al. 1992) Even for systems with a number of triple-value zones, the first guess 
of final and initial particle position pairings produces rapid convergence. One can then 
permute sf into sf , resulting in a well-motivated set of boundary constraints. 



2.2. Computing Trajectories in Redshift Space 

We will now consider the simultaneous determination of the orbits of interacting 
particles when the boundary constraints are given in redshift space. Let us say that we 
have an isolated, uniform grid of particles at a = 0, with positions given by {q^}, and at 
t = t , those particles are "observed" at redshift coordinates, {sf}. We write down the 
trajectories of the particles as the sum of a part given by linear perturbation theory, and 
coefficients times basis functions. However, unlike the discussion of PLA in real space (GS), 
redshift space gives us heterogeneous constraints on our basis functions. In real space, the 
basis functions were constrained such that: 

/„(*„) = ; lima 2 /n(t)=0 (21) 

a— >0 

In GS, we showed that these constraints can be satisfied by using 

f n = E b nm D(t) m [D - D(t)\ , (22) 

m=n 

where b lm = 5^ m , and the higher order coefficients are based on fitting to an N-body 
simulation. 

In redshift space, however, we need to define a slightly different set of basis functions, 
f n (t). In this case the basis functions along the line of sight must satisfy 

/„(t ) + ^M = ; lima 2 / n (t) = 0, (23) 
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in order that the varying the coefficients of f n do not change the corresponding radial 
redshift space coordinate. 

To satisfy these constraints, we introduce a complementary set of basis functions to 
those introduced in real space: 

Ut) = f n {t) - . (24) 

ao 1 + n 

It can be shown that these functions satisfy equation (p3|) if the real space basis functions 
described above are used. As with the unaccented basis functions, only the first function 
goes linearly or slower at early times. Using these basis functions, we are able to describe 
the trajectory of any particle as: 

x «(t)=xf ,)a + D(t)(sf-xf ,)a )+J2CUn(t) ; a = 1,2 (25) 

n 

= + n ; j ;,o ~ c x ; f I + E <%/*(*) ; « = 3 , w 

where xf^ a are the set of some physically self-consistent orbits, as output from an N-body 
code. Of course, we may also set xf^ a = qi. 

In GS, we showed that by simultaneously minimizing the action, S, for each coefficient, 
Cf n , then the equations of motion of the particles are necessarily satisfied, and hence, we 
may find these orbits. That is, 



OS fto 



f °dt - *£) + 

Jo 1 



4(a 2 xp) + ^ 0) - ()<> > 



dt K 1 ' dx? dxf 



(27) 



We will solve equation ( |2"7| ) for all possible basis functions simultaneously, by 
determining the coefficients such that the kernel 

g,( t ) = -*&p. + (V# - V*) (28) 

vanishes at all times, and for all particles. 

In order to do this, we can use the metric: 

X 2 = Y,f dt \W(tMt)\\ (29) 

i 

where W(t) is an arbitrary weighting function. By minimizing X 2 , we find the set of 
trajectories which come closest to satisfying the equations of motion. 
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Perturbations of the basis functions may be approximated by: 



5C 7 ~ 



. ' dtW(t) 2 gf old \t)A^(t) 

m-y,nf3 



n(3 



(30) 



where 

A „ , rt _*>»?(*) s d\a'xf ) ut)\ ay, 

A„„,,(0 = g^- = -** m ^« 4 > ■ < 31 > 

Where (3 = 3, the basis functions, /„(£) should be replaced by /„(£). 

We illustrate the results of this method on the Zel'dovich pancake in Figure [I] and 
Figure For this, we have used 4 basis functions and two iterations. For each iteration, 
we took the constraint pairs, {q^, sf }, and used PLA to compute the best fit full trajectory. 
We then evaluated the positions and velocities of the particles at a ~ 0.01, and ran the 
particles through the PM code again. 

The long dashed lines in panels c) and d) of Figure [I] represent the density field of 
the second iteration in redshift and real space, respectively. While by the nature of the 
constraints we would necessarily expect the redshift space density field to converge to 
the "true" field, we have no such guarantee in real space. Nevertheless, the real space 
density field does seem to give a somewhat better fit than the initial estimate given by the 
VELMOD-like approach, especially around the edge of the peak, PLA gives a smoother 
edge. 

Additionally, even though the peak is nonlinear in the sense that 5 > 1 in real space, 
PLA is able to very successfully generate an initial density field. Figure [l|a) and b) show the 
redshift and real space density field determined by PLA, as well as the true initial density 
field. PLA is able to determine the amplitude of the initial peak to about 10%. However, 
it should be noted that PLA may erroneously generate too much small scale power. Given 
some a priori knowledge of the power spectrum, however, one may use a power-preserving 
filter like the one described in GS in order to appropriately smooth the field. In this 1-d 
case, in which we only expect a single mode, using such a filter would be gratuitous. 

Finally, in Figure we show the relationship between real and redshift space 
coordinates as determined from the output of the PLA code. Note that the PLA points turn 
over more smoothly than our initial guess points. This is due to the fact that PLA assumes 
the field to be evolving from an initially uniform particle field, while the VELMOD-like 
approach makes no such assumption. 
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3. A High Resolution Test of the Code 

To illustrate PLA's success as a reconstruction scheme somewhat into the nonlinear 
regime, we have run a high resolution simulation, with N p = 128 2 , N g = 256 3 , and with 
L = 800/i _1 Mpc. The cosmology used is Qm = 0.3, Q\ = 0.7, and as = 0.7. We then apply 
PLA under the assumption that the cosmology was known to reconstruct the field. 

In Figure |3| we show a density contour of 5 = 0.7 for the observed redshift space density 
field, smoothed with a Gaussian filter with a radius of Mpc, and a similar plot, but for 
the field resulting from our PLA reconstruction. Recall that the underlying particle field 
for the latter is based on running the reconstructed initial conditions through a PM code 
and taking the smoothed density field. A visual inspection demonstrates that the two fields 
are virtually identical. 

We plot a similar comparison in Figure f|, except with a smoothing radius of only 4/i _1 
Mpc, and a density contour of 5 = 2.1. On these scales, too, the fields are very similar. 
The reconstructed field, however, seems to differ somewhat on small scales. If we had an a 
priori model of the power spectrum, this small scale power could be suppressed numerically. 

More quantitatively, Figure |5] shows the fit between the reconstructed and true field as 
a function of scale. Narayanan & Croft (1999) provide the goodness of fit metric: 

{ ' j "£[5i(k) 2 + 5 2 (k)p • { } 

In this case, #i(k) and ^(k) are the Fourier transforms of the true and reconstructed 
redshift density fields. For perfect matching on a particular scale, this metric goes to zero. 
For uncorrelated fields, it goes to one. 

As Figure |5| illustrates, the fit for all four comparisons is very good even into the 
nonlinear regime (k ~ OAh Mpc- 1 ). The best fit was found for the final conditions in 
redshift space, since this was the actual set of observations to be matched. However, 
it is shown that the real space field at z=0 is also reconstructed quite well, as are the 
corresponding initial conditions. In particular, the relevant scale is that for which A 2 = 0.5 
for each comparison. In this case, the real and redshift initial conditions are well matched 
down to a scale of 18.6 and 19.5/i -1 Mpc, respectively, and the real and redshift space final 
conditions are matched down to scales of 15.1 and 12.7ft. -1 Mpc, respectively. 



Another measure of the quality of the reconstruction is the isotropy of the reconstructed 
real space density field. We may do this by decomposing the power spectrum of the field 
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into Legendre Polynomials (e.g. Hamilton, 1992): 

p(k) = e nwm) , (33) 

l even 

where /i is the cosine of the angle between the direction vector, k, and the line of sight, 
Vi(fi) are the Legendre polynomials, Pi(k) are the azimuthally-averaged multipole expansion 
of the power spectrum. Taking the inverse Legendre transform, we find (Cole, Fisher and 
Weinberg 1995): 

2/ + 1 r 1 r 2n 
Pi(k) = \ dfi / d<t>P{k)Vi{fj,) . (34) 

Air J-i Jo 

Here, we will only be using the monopole and quadrupole moments, which, as a reminder 
are Vo(fi) = 1 and ^(/-O = (3yU 2 — l)/2. We may thus measure the isotropy of the 
distribution by computing the ratio: 

<?.(*) s m . <35) 

Since the real universe is isotropic, an accurately computed reconstruction should have a 
quadrupole equal to zero on all scales. 

Figure |] shows the quadrupole ratio as a function of scale for both 1 and 2 iterations. 
Two things are clear from this plot, however. First, the overall isotropy does improve with 
additional iterations. Secondly, there is a systematic effect in generating these anisotropies 
which is almost certainly caused by an anisotropic noise term. In the appendix, we use the 
quadrupole and hexadecipole moments to illustrate that this effect is dominated by noise, 
rather than by a systematic underestimate of the radial velocity term, for example. 

This noise term comes out of the reconstruction scheme, itself. In the previous section, 
we discussed how one goes about approximately rewinding the trajectory of a particle in 
order to determine its initial constraint. However, there was an assumption that linearity 
approximately held. As structure gets more and more nonlinear, this assumption will 
fail to hold, and the particle matching technique will break down. Future advances in 
reconstruction methods will have to take this into account. It may be possible to apply 
PLA in an iterative and statistical way in order better do this matching. 



4. Applications: Breaking the Bias Degeneracy 
4.1. Motivation 



Though the reconstruction of nonlinear fields is interesting in its own right, we have 
begin this investigation into redshift space for the purpose of finding out something about 
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cosmology. As a motivation for this sort of reconstruction analysis of redshift surveys, 
we address the M -6 degeneracy in redshift surveys, and discuss how the degeneracy 
might be broken without recourse to outside dynamical estimates of Qm- We will test the 
effectiveness of using PLA to break the degeneracy. 

We begin by introducing the problem as it appears in the linear regime. Excellent 
recent reviews of this topic is given by Hamilton (1998) and by Strauss & Willick (1995), 
and we will therefore present only an overview of the linear biasing problem. 

Let us first suppose that we have the full velocity field information about a group 
of almost uniformly distributed particles. Let us further suppose that these particles are 
biased with respect to some true underlying CDM field, such that: 

5 g (x,t ) = b5(x,t ) , (36) 

where the unsubscripted 5 is the CDM density field, and 5 g (for galaxies) represents the 
density field of some biased tracer of the mass. Note that a straight linear biasing model 
may be replaced with any deterministic, local, and monotonically increasing function of 
5(x, to) without changing the essence of the discussion or the biasing problem in general. 
Blanton (1999) provides an excellent review of various types of biasing models. 

With linear bias greater than unity the velocity field is of lower amplitude than one 
would directly infer from a measurement of 5 g . From (^), equation (||) may be replaced by: 

V-v(x,t ) = -P5 g (x,to) , (37) 

where 

b b 



and 

v = f(n M ,n A ,t ) P . (39) 

Thus, the divergence of the velocity field and density field are related by the same 
proportionality constants for all combinations of Qm and b which yield the same (3. This is 
the crux of the degeneracy problem. 

In a redshift survey, we do not actually know the divergence of the velocity field, but 
rather must infer it through anisotropics in the redshift density field. Kaiser (1987) shows 
that in the linear regime, there is a straightforward relationship between the real and 
redshift space density fields in k-space in the d.o.a.: 



^(k) = s^ = ^(k)[i + /V] , 



(40) 
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where S is the linear redshift distortion operator, <5*(k) is the Fourier transform of the 
galaxy density field in redshift space, and 5 x (k) is the Fourier transform of the underlying 
CDM density field in real space. 

From this relationship, a redshift space power spectrum field may be compared to the 
real space power spectrum: 

P5(k)^(k)| 2 = P(k)[l + /V? • (41) 



We can further decompose the power spectrum field into Legendre polynomials (equation 
Hamilton 1992). We then take the inverse Legendre transform (equation |34|; Cole, Fisher 
and Weinberg 1995) and compute the quadrupole moment. 



Taking the Legendre expansion of the form of the linear density field (equation |33|) 
above, we find that the moments of the biased, redshift space distribution may be related 
to the underlying real space distribution as follows (Hamilton 1998): 



P°,o( k ) = [l + \(3+\(3 2 )p{k) (42) 
P:jk) = [\ ? • I ?) />(/■•} • (43) 



3' 7 



Though the underlying power spectrum is not known, this quadrupole ratio may be 
computed as: 

W = SiW = (44) 

Since the multipole expansion may be computed directly from the observed redshift space 
density field, it is clear that in the linear regime, the quadrupole ratio is degenerate for a 
particular value of (3. 

In addition, (3 may be estimated through the use of distance-velocity comparisons. In 
practice, measurements of (3 are still quite difficult due to the noisy data involved. Willick 
(2000) gives a review of current estimates, and finds a value of (3 = 0.5 ± 0.04 from the 
IRAS velocity field and the Mark III Catalog used to estimate distances. This is relatively 
unchanged from the earlier estimate based on measured redshift-space anisotropy by Cole, 
Fisher, & Weinberg (1995) of (3 = 0.54 ± 0.3. Ballinger et al. (2000) estimate /3 = 0.4 ± 0.1 
for the IRAS 0.6 Jy PSCz survey from analysis of anisotropies. Similar values are found 
from the Optical Redshift Survey catalog (Baker et al. 1998). 

After estimating f3, how does one estimate Qm without recourse to direct mass 
estimates? We answer this by first pointing out that a correctly reconstructed real space 
density field will be perfectly isotropic. That is, Q x (k) = on all scales. Nonlinearities 
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is the evolution of the density field will mean that the linear redshift distortion operator 
will no longer be valid on all scales. By reconstructing fields using PLA for different 
combinations of Qm and 6, and measuring the quadrupole moment for the reconstructed 
field, we can determine the "true" cosmology as that which minimizes the anisotropy. 



4.2. Simulations 

4-2.1. Constraining Bias 

In order to test this approach, we have run two simulations, each extending only into 
the mildly nonlinear regime. We have found that simulations which contain highly nonlinear 
structure do not effectively differentiate between different cosmologies due to the excessive 
noise and difficulty of doing the particle matching on small scales. 

The two simulations were each run with j3 = 0.5 and a boxsize of L = 1000/i _1 Mpc, 
the first with Qm = 0.3, Q\ = 0.7, and 6 = 1, and the second with Qm = 1, = 0, 
and 6 = 2. Each simulation was run with 64 3 particles, and 128 3 gridcells, and the PLA 
reconstruction used 4 basis functions. The "observations" of these simulations were the 
redshift-space density fields in the d.o.a. After reconstructing the initial density field for a 
particular assumed bias, we ran the initial conditions through a PM code, and measured 
the quadrupole moment in real-space, where it should vanish. The "mean" quadrupole 
moment is estimated as: 

/ n 2\ _ I dk Ql(k) k 2 e k l kl 

Wx) ~ fdkk*e-k/*i ' 1 j 

where ki is a limiting scale at which point grid and/or nonlinear effects will become 
important. We have assumed ki = 20 x 2ir/Li )OX , but found similar results for ki = oo. The 
peak of this contribution occurs around 3ki, or on a physical scale of ~ lQh~ l Mpc, the 
nonlinear scale. The results of each model tested are shown in Figure for the Qm = 0.3 
simulation, and Figure [8] for the Qm = 1 simulation. 

In each case, we find that the best fit value of Qm corresponds to the actual value of 
Qm used in the simulation. That is, by tracing the detailed evolution of a mildly nonlinear 
field, PLA can effectively break the j3 degeneracy. 

In estimating this effect from a real survey, we would have to Monte Carlo observations 
based on the survey geometry, selection function, and the like, in order to estimate Q m and 
its errors. While direct error estimation is difficult with only two realizations, the results 
are quite suggestive that this will be an effective way to constrain Q M directly. Susperregi 
(2000) uses an Eulerian least action code to similarly show that the bias degeneracy may 
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be broken through accurate reconstruction. 

As a final test of PLA as a redshift space reconstruction scheme, we show that one 
may obtain somewhat better estimates of (3, itself, from the assumed isotropy of the 
reconstructed real space density field than from the redshift space anisotropy. This is 
illustrated in Figure U, in which we show a comparison between Q x {k) of the reconstructed 
field for simulation 1, assuming the correct value of Qm — 0.3, and the corresponding 
quadrupole moment residuals, Q s {k) — of the observed redshift space field. Note that 
the term is simply that which one would estimate from an assumption of the correct 
value of p. 

On large scales, these two statistics are almost identical. However, on smaller 
scales, when nonlinearities begin to become important, the two statistics both exhibit 
anisotropies. The redshift space field, however, becomes anisotropic on larger scales. Since 
the corresponding uncertainties in (3 are approximately proportional to one over the square 
root of the number of modes probed, we may relate the expected uncertainties from the 
reconstructed field to that from the redshift field as: 

„PLA n PLA d(i (j PLA\-3/2 
°P Q 8Q Z v^maxJ 

where the superscript "Z" refers to the estimate from the redshift space field and the 
superscript "PLA" refers to the estimate from the reconstructed real space field. For 
Qm = 0.3, all of the partial derivatives are almost exactly one. Moreover, estimates of the 
scatter in the quadrupole estimates show that &p LA — 0.9a^. Finally, since Figure |9| shows 
^max 4 — l-25fc max , the approximate relation between the uncertainty in the bias between the 
two errors is <Jp LA — 0.65crf . Thus, we are able to generate a somewhat better constrained 
estimate of j3 from the reconstructed density field, as compared to the observed field. 



5. Future Prospects 

This paper has discussed the problem of reconstructing the underlying real space CDM 
density field and its evolution from galaxy redshift surveys under rather idealized conditions. 
We have assumed full sampling, a constant linear bias relation with no morphological 
segregation, the distant observer approximation, no errors in measurement, and a very 
regular geometry. We have shown that theoretically even a mildly nonlinear field can be 
reconstructed using PLA to break the bias degeneracy. In actual observations almost none 
of these assumptions will hold. We would like to end this work with a brief discussion of 
how these effects might be appropriately modeled, and mention a few possible candidates 
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of real surveys to which the PLA method might be applied. 

Throughout, we have assumed a cubic geometry. In simulating a realistic survey, a 
mask must be applied such that statistics may be correctly computed for the true survey 
volume. Additionally, for many observational samples of interest, the distant observer 
approximation no longer describes the system adequately. Over the course of its lifetime, 
a galaxy may have traversed a significant angle in the sky, making the use of Cartesian 
coordinates difficult. An obvious solution to this problem is to write the PLA equations in 
spherical coordinates. 

An additional concern in the application of PLA to realistic observations is that the 
biasing model that we have used here is a strictly linear one, and observational evidence 
suggest that biasing may be much more complex (e.g. Blanton 1999 and references therein). 
In fact, we have only used linear bias in this work for its simplicity. PLA would be 
applicable to any model of deterministic bias, even one which had an explicit model of 
morphological segregation. It is not clear how one might effectively reconstruct a field under 
the assumption of a significantly stochastic bias model. 

Even under the simplest of approximations, realistic surveys are not cubic, not 
necessarily contiguous, and generally have nonuniform selection functions. The issues of 
contiguity and geometry are related, in that in both we need to approximate the density 
field outside the survey volume in order to correctly estimate the potential field within. 

Lahav et al. (1994) apply one such technique to the IRAS 1.2 Jy survey (Strauss et 
al. 1992; Fisher et al. 1995) in which expansion of the density field in spherical harmonics 
is used to reconstruct the field outside the survey volume. In essence, this is very similar 
to assuming a particular autocorrelation function, and generating an outside field based 
on a truncated form of that function and the observed field near the edges. In many 
respects, this is quite similar to the sort of reconstruction done using the "Constrained 
Initial Conditions" technique by Hoffman & Ribak (1991,1992), since both use observed 
the observed autocorrelation function to build realistic external fields around observed 
structure. 

The final issue in real surveys concerns non-uniformity and noise within the survey 
volume. However, for flux limited surveys such as the IRAS 1.2 Jy and PSCz (Saunder 
et al. 2000) surveys, at large distances, shot noise begins to dominate calculations of the 
density field. Since the observed density field is derived from an incomplete sampling of a 
finite number of discrete points, the uncertainties in the observed density field behaves like 
a Poisson statistic, with a oc 1/ ^ N ga i . In order to correctly anticipate the effects of shot 
noise when running simulations of a particular survey, a random component needs to be 
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added to the observed density field. 

Moreover, a non-uniform selection function, coupled with observations in redshift space 
results in Malmquist bias. That is, the "true" selection function is based on observed fluxes, 
and hence is a function in real space. The observations, however, are of densities in redshift 
space. Strauss & Willick (1995) give an excellent review of how these issues are dealt with 
in reconstructing density fields. 

These considerations are all with an eye toward learning about cosmology from redshift 
surveys. For example, the IRAS 0.6 Jy PSCz Survey (Saunders et al. 2000) is an especially 
promising recent candidate for analysis, as it is publicly available and has nearly full sky 
coverage. Hamilton, Tegmark, & Padmanabhan (2000) have already estimated (5 = 0.411^2 
from this survey, and Ballinger et al. (2000) find a similar result of /3 — 0.4 ±0.1. Both 
methods used only linear theory, however. Valentine, Saunders, & Taylor (2000) do a 
somewhat higher order reconstruction, by using the PIZA method, and find a best fit to 
the survey with a slightly higher value of f3 ~ 0.5. In testing our code, we have found 
that one can get a better measure of (5 by using PLA to reconstruct a density field under 
some fiducial cosmology, and comparing cosmologies to see which produce the minimum 
anisotropy in the real space field. We have showed that ideally, PLA can be used to 
discriminate between "degenerate" pairs of bias and Qm- Finally, we showed that PLA 
produces a means by which uncertainties in the measurement of f3 itself can be reduced. 

Another interesting prospect is the application to PLA to large distance redshift 
surveys, since one of the byproducts of PLA is the real space density field. Nusser et al. 
(2000) use the D n — a relation of the ENEAR redshift-distance survey (da Costa et al. 
2000) as test particles within the PSCz survey, much as we would wish to do using PLA. A 
reconstruction was done using the method described by Nusser & Davis (1995). Predicted 
distances from the reconstructed field can then be compared with the estimated distances 
from the ENEAR survey. They also estimated (3 ~ 0.5. 

While the PSCz survey contains ~ 15, 000 redshifts, the current generation of redshift 
surveys is producing an even greater opportunity to measure statistical and global properties 
of the universe. When complete, the SDSS redshift survey (York et al. 2000), will produce 
~ 10 6 galaxy redshifts, and will cover a quarter of the sky out to Petrosian magnitude, 
r' p = 17.7 toz~ 0.15. The 2dF survey (Colless 1999) will ultimately measure redshifts over 
a quarter of a million galaxies, out to bj = 19.5. 

Reconstruction of fields from these enormous datasets will prove a significant 
computational challenge. However, it is well worth it, as PLA can yield insight into the 
underlying power spectrum, bias, and cosmological parameters. 
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A. The Effects of Noise on M ultipole Moments 



In §3, we found that the reconstructed real space density field of the high-resolution 
simulation had a significant and seemingly systematic quadrupole moment at small scales. 
The question which arises out of this is, does this quadrupole represent a systematic error in 
the reconstruction along the line of sight (such as would occur, for example, if the assumed 
value of f3 were incorrect), or does it represent a stochastic term? 

To examine this question, let us consider the following form of a reconstructed field 
with a very simple error term: 



<5(k) = 5(k) 1 + fi 2 B(k) + fi 2 e(k) 



(Al) 



where k is a Fourier space component, 5(k) is the true real space field we are attempting to 
reconstruct, <5(k) is the reconstructed form of the field, fi 2 is the cosine of the angle between 
the k and the line of sight, B is a systematic, "bias" term, and e(k 2 ) is an anisotropic 
random error drawn from a N (0 , a 2 (k)) distribution. 

If the corresponding real space density field of redshift space observations have been 
perfectly reconstructed and contain no noise, the reconstructed field is simply equal to 
the true underlying field, and thus will be perfectly isotropic in A;-space. However, let us 
imagine that a field contains no noise, but the reconstruction is such that we have assumed 
that 5(k) = 5 s (k), or the real space density field is the same as the redshift space field. 
Under those circumstances, B(k) = ft, hence our terminology. 

Finally, let us consider a more general case, one in which we wish to test for a 
systematic form of B(k) and for the existence of a random noise component. In that case, 
the reconstructed three-dimensional power spectrum may be written as: 

/P(k)N 



i + 2/2 2 B(ky + ^B(ky + /iV(k) 



(A2) 



If we then decompose these terms into multipole moments (see §5.4.1) and assume that 
the noise and systematic terms are simply scale dependent, we find: 



Po(fc) = P(k) 
P 2 (k) = P(k) 
P 4 (k) = P(k) 



1 + 2 -B{k) + \B{k) 2 + \<j 2 {k) 
3 5 5 

js(fc) + ^B(k) 2 + A -a\k) 



^B(k) 2 + ^(*) 



(A3) 
(A4) 
(A5) 



We may then look at the relationship between the quadrupole ratio, Q(k) = P 2 (k) / P (k) 
and the hexadecipole ratio, H(k) = P 4 (k)/P (k). We have two extremes. In the case of no 
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anistropic noise component, each ratio is simply a parametric function of B(k), and thus, a 
straightforward relation between the two may be plotted. 

If, on the other hand, the anisotropic noise term dominates, we find the relation: 

H NOISE (k) = 2 -Q NOISE (h) . (A6) 





We illustrate this in Figure [TC]. We show that for an observed redshift density field, the 
deterministic bias term dominates. Though we do not necessarily have the correct form of 
the anisotropic error term, this error which would seem to be at the root of the systematic 
small scale quadrupole in the high-resolution simulation. 

Since we may take as a prior that the universe is inherently isotropic, regularization or 
iterative techniques might be employed in future reconstruction schemes which explicitly 
find a fully isotropic real space solution. 




Fig. 1. — A simulation of a quasi-linear Zel'dovich pancake. In each panel, the solid 
lines show the "true" real- and redshift space density fields before and after running the 
corresponding particles through an N-body code. Panels a) and b) show the initial density 
fields in redshift and real space, respectively. The long-dashed lines represent the initial 
density field computed by using PLA. Panel c) shows the evolved field in real space. The 
solid line shows the true field, while the dotted line shows the first estimate of the field from 
the redshift space distribution, and the short dashed line shows the converged value. The 
long dashed line shows the density field computed using PLA. Panel d) shows the evolved 
field in redshift space. The dotted line denotes the density field corresponding to the redshift 
positions generated by using the laminar flow density matching method described in the text. 
The long dashed line shows the density field computed by using PLA. 
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Fig. 2. — The real space-redshift space coordinate relationship for a Zel'dovich pancake. 
The solid line shows the relationship as given by the output of an N-body simulation. The 
squares show the relation as given by the VELMOD-like scheme used to approximate the 
real-space density field. Finally, the solid points show the output from PLA. 
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Fig. 3. — A comparison of the smoothed redshift space density fields of the output of 
the resolution simulation described in the text, and the reconstructed field based on those 
observations. A contour of 5 = 0.7, smoothed with a Gaussian filter of r = %h~ x Mpc is 
shown. Note the high level of agreement between structures in the true and reconstructed 
fields. 
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Fig. 4. — As in the previous figure, but with a Gaussian smoothing radius of 4/i _1 Mpc, and 
a density contour of 8 — 2.1. 
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Fig. 5. — The Fourier difference statistic for the high resolution simulation. This compares 
the true initial and final conditions in real and redshift space to those generated by using the 
reconstruction scheme using only the observations of the final redshift space density field. 
The dotted line shows the first iteration of PLA, and the solid line, an additional iteration. 
Note that a physical scale of i6/i -1 Mpc (the nonlinear scale) occurs at log 10 (/c) = —0.4. 
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Fig. 6. — The quadrupole moment ratio as a function of \k\ for the reconstructed high- 
resolution real space density field simulation. Solid squares represent the results from one 
iteration of PLA, while solid circles are the results from a second iteration. 
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Fig. 7. — The weighted-average real space quadrupole ratio for the reconstructed density 
fields with various assumed cosmologies in the VL M = 0.3 simulation. Note that the 
quadrupole ratio is minimized for Qm = 0.3, indicating that PLA reconstruction effectively 
breaks the f3 degeneracy. 




Fig. 8. — As in the previous figure, but for the Qm — 1 simulation. Here, the quadrupole 
ratio is minimized for Qm = 1, again indicating the reconstruction method has broken the 
degeneracy. 
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Fig. 9. — A comparison of the scale dependent quadrupole moment ratio for the reconstructed 
real space density (filled circles), and the "observed" redshift space density (open squares), 
with the linear quadrupole term subtracted. Each estimate was done using the Qm = 0.3 
low resolution simulation. Note that while both statistics are approximately zero at large 
scales, there is a systematic divergence at small scales. This turnoff occurs at larger scales 
for the direct redshift estimate than for the estimate based on the reconstructed field. 
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Fig. 10. — Scatter plots relating the quadrupole and hexadicpole ratios in the high-resolution 
simulations discussion in §3. The solid line shows the expected relation for a systematic 
anisotropy (the /3-like term), while the dashed line shows the expected relation for our 
simplified anisotropic noise model. Panel a) shows the ratios as observed in the redshift space 
density field. Note that the systematic effect dominates, since the anisotropics are completely 
due to redshift space distortions. Panel b) shows the ratios for the reconstructed field. The 
hexadecipole ratio is much higher than what would be expected from a deterministic effect. 





